Novel basin modelling concept for simulating deformation from mechanical compaction using level sets

As sedimentation progresses in the formation and evolution of a depositional geologic basin, the rock strata are subject to various stresses. With increasing lithostatic pressure, compressional forces act to compact the porous rock matrix, leading to overpressure buildup, changes in the fluid pore pressure and fluid flow. In the context of petroleum systems modelling, the present study concerns the geometry changes that a compacting basin experiences subject to deposition. The purpose is to track the positions of the rock layer interfaces as compaction occurs. To handle the challenge of potentially large geometry deformations, a new modelling concept is proposed that couples the pore pressure equation with a level set method to determine the movement of lithostratigraphic interfaces. The level set method propagates an interface according to a prescribed speed. The coupling term for the pore pressure and level-set equations consists of this speed function, which is dependent on the compaction law. The two primary features of this approach are the simplicity of the grid and the flexibility of the speed function. A first evaluation of the model concept is presented based on an implementation for one spatial dimension accounting for vertical effective stress. Isothermal conditions with a constant fluid density and viscosity were assumed. The accuracy of the implemented numerical solution for the case of a single stratigraphic unit with a linear compaction law was compared to the available analytical solution [38]. The multi-layer setup and the nonlinear case were tested for plausibility.     

A comparative simulation study of coupled THM processes and their effect on fractured rock permeability around nuclear waste repositories

This paper presents an international, multiple-code, simulation study of coupled thermal, hydrological, and mechanical (THM) processes and their effect on permeability and fluid flow in fractured rock around heated underground nuclear waste emplacement drifts. Simulations were conducted considering two types of repository settings (1) open emplacement drifts in relatively shallow unsaturated volcanic rock, and (2) backfilled emplacement drifts in deeper saturated crystalline rock. The results showed that for the two assumed repository settings, the dominant mechanism of changes in rock permeability was thermal–mechanically induced closure (reduced aperture) of vertical fractures, caused by thermal stress resulting from repository-wide heating of the rock mass. The magnitude of thermal–mechanically induced changes in permeability was more substantial in the case of an emplacement drift located in a relatively shallow, low-stress environment where the rock is more compliant, allowing more substantial fracture closure during thermal stressing. However, in both of the assumed repository settings in this study, the thermal–mechanically induced changes in permeability caused relatively small changes in the flow field, with most changes occurring in the vicinity of the emplacement drifts.     

Hybrid analytical and finite element numerical modeling of mass and heat transport in fractured rocks with matrix diffusion

Quantification of mass and heat transport in fractured porous rocks is important to areas such as contaminant transport, storage and release in fractured rock aquifers, the migration and sorption of radioactive nuclides from waste depositories, and the characterization of engineered heat exchangers in the context of enhanced geothermal systems. The large difference between flow and transport characteristics in fractures and in the surrounding matrix rock means models of such systems are forced to make a number of simplifications. Analytical approaches assume a homogeneous system, numerical approaches address the scale at which a process is operating, but may lose individual important processes due to averaging considerations. Numerical stability criteria limit the contrasts possible in defining material properties. Here, a hybrid analytical–numerical method for transport modeling in fractured media is presented. This method combines a numerical model for flow and transport in a heterogeneous fracture and an analytical solution for matrix diffusion. By linking the two types of model, the advantages of both methods can be combined. The methodology as well as the mathematical background are developed, verified for simple geometries, and applied to fractures representing experimental field conditions in the Grimsel rock laboratory.     

ESR Dating of the Evolution of the Shuanghu Basin in the Northern Tibetan Plateau

The Shuanghu basin is a NE-trending rift basin bounded by NE-striking normal faults and NW-striking shear-extensional faults of the northern Tibetan Plateau. Four samples from calcite veins in marginal faults and one sample from mudstone (S-3) were collected for dating the evolution of the Shuanghu basin by using the ESR spectrograph of EXM-type. Ages were calculated according to the close-equilibrium model on the basis of the measured ESR signal spectra of samples, providing good chronological information. It is known from the ESR dating that the extensional faulting and rifting of the Shuanghu area began at 4.92 Ma B.P., followed by regional folding in 3.56-1.36 Ma, NW-striking faulting in 0.60 Ma and normal faulting in 0.024 Ma in the Shuanghu basin.     

Insect communities of Phragmites habitats used for sewage purification: Effects of age and area of habitats on species richness and herbivore-parasitoid interactions

In 1994, 28 sewage purification plants based on common reed Phragmites australis were analyzed, to study fragmentation effects on diversity and on food web structure of insect communities. Purification plants differed in age (2 to 11 years) and area (10 to 2,500 m²). Both age and area were significantly correlated with species richness of herbivores. Populations of the two most abundant herbivores, the gall makers Lipara pullitarsis (Diptera, Chloropidae) and Giraudiella inclusa (Diptera, Cecidomyiidae), increased significantly with age and size of the Phragmites habitats, i.e. the purification plants. Incidence curves showed that only habitats older than 6 years and larger than 100 m² have a 50% probability of finding these two most abundant gall makers, and can correspondingly be considered to contribute to insect species conservation. Natural enemies were even more susceptible to habitat fragmentation than their phytophagous hosts or prey, thereby hampering possible biocontrol. Percent parasitism of both the Lipara and Giraudiella gall makers significantly increased with habitat age.     

Is the North Atlantic Oscillation reflected in Scandinavian glacier mass balance records?

The North Atlantic Oscillation (NAO) is one of the modes of climate variability in the North Atlantic region. The atmospheric circulation during the winter season in this region commonly displays a strong meridional (north–south) pressure contrast, with low air pressure (cyclone) centred close to Iceland and high air pressure (anticyclone) near the Azores. This pressure gradient drives the mean surface winds and the mid‐latitude winter storms from west to east across the North Atlantic, bringing mild moist air to northwest Europe. The NAO index is based on the difference of normalised sea‐level pressures (SLP) between Ponta Delgada, Azores and Stykkisholmur, Iceland. The SLP anomalies at these stations are normalised by division of each monthly pressure by the long‐term (1865–1984) standard deviation. Interannual atmospheric climate variability in northwest Europe, especially over Great Britain and western Scandinavia has, during the last decades, been attributed mainly to the NAO, causing variations in the winter weather over the northeast North Atlantic and the adjacent land areas. A comparison between the NAO index and the winter (December–March) precipitation between ad 1864 and 1995 in western Norway shows that these are strongly linked (correlation coefficient 0.77). Variations in the NAO index are also reflected in the mass balance records of glaciers in western Scandinavia. The NAO index is best correlated with mass balance data from maritime glaciers in southern Norway (e.g. Ålfotbreen R² = 0.51). The record of Holocene (last ca. 11 500 cal. yr) glacier variations of maritime glaciers in western Scandinavia is thus a proxy of pre‐instrumental NAO variations. Copyright © 2000 John Wiley & Sons, Ltd.     

Gamma rays from colliding winds of massive stars

Colliding winds of massive binaries have long been considered as potential sites of non-thermal high-energy photon production. This is motivated by the detection of non-thermal spectra in the radio band, as well as by correlation studies of yet unidentified EGRET γ-ray sources with source populations appearing in star formation regions. This work re-considers the basic radiative processes and its properties that lead to high energy photon production in long-period massive star systems. We show that Klein–Nishina effects as well as the anisotropic nature of the inverse Compton scattering, the dominating leptonic emission process, likely yield spectral and variability signatures in the γ-ray domain at or above the sensitivity of current or upcoming gamma ray instruments like GLAST-LAT. In addition to all relevant radiative losses, we include propagation (such as convection in the stellar wind) as well as photon absorption effects, which a priori can not be neglected. The calculations are applied to WR 140 and WR 147, and predictions for their detectability in the γ-ray regime are provided. Physically similar specimen of their kind like WR 146, WR 137, WR 138, WR 112 and WR 125 may be regarded as candidate sources at GeV energies for near-future γ-ray experiments. Finally, we discuss several aspects relevant for eventually identifying this source class as a γ-ray emitting population. Thereby we utilize our findings on the expected radiative behavior of typical colliding wind binaries in the γ-ray regime as well as its expected spatial distribution on the γ-ray sky.     

Use of a parsimonious rainfall–run‐off model for predicting hydrological response in ungauged basins

A new parameter parsimonious rainfall–run‐off model, the Distance Distribution Dynamics (DDD) model, is used to simulate hydrological time series at ungauged sites in the Lygne basin in Norway. The model parameters were estimated as functions of catchment characteristics determined by geographical information system. The multiple regression equations relating catchment characteristics and model parameters were trained from 84 calibrated catchments located all over Norway, and all model parameters showed significant correlations with catchment characteristics. The significant correlation coefficients (with p‐value < 0.05) ranged from 0.22 to 0.55. The suitability of DDD for predictions in ungauged basins was tested for 17 catchments not used to estimate the multiple regression equations. For ten of the 17 catchments, deviations in Nash–Sutcliffe efficiency (NSE) criteria between the calibrated and regionalised model were less than 0.1, and for two calibrated catchments within the Lygne basin, the deviations were less than 0.08. The median NSE for the regionalized DDD for the 17 catchments for two time series was 0.66 and 0.72. Deviations in NSE between calibrated and regionalised models are well explained by the deviations between calibrated and regressed parameters describing spatial snow distribution and snowmelt respectively. The quality of the simulated run‐off series for the ungauged sites in the Lygne basin was assessed by comparing flow indices describing high, medium and low flow estimated from observed run‐off at the 17 catchments and for the simulated run‐off series. The indices estimated for the simulated series were generally well within the ranges defined by the 17 observed series. Copyright © 2014 John Wiley & Sons, Ltd.     

Multipole borehole acoustic theory: Source imbalances and the effects of an elastic logging tool

In recent years the emphasis in acoustic logging has been shifting from the wireline to the Logging While Drilling (LWD) environment, the latter being far different from the former in that both tool rigidity and tool radius are considerably greater. In this paper we present a generic mathematical formulation for the multipole borehole acoustic measurement (alternate and equal polarity case), including a detailed analysis on the effects of multipole source amplitude imbalances. It is shown that source imbalance induced mode contaminants have excitation amplitudes that are scaled by the sum of the relative source imbalances between diametrically opposed sources. Furthermore it is shown that for mode contaminants with odd modal number there can be a significant offset in the associated directivity pattern, even at low levels of source imbalance. However, it is also shown that source imbalance induced mode contaminants can be completely eliminated if the multipole source is accompanied with a ‘vertically’ (but not azimuthally) offset multipole receiver. Mathematically, it is demonstrated (for a centered tool) that a polarity weighted stack of these multipole receivers completely eliminates the source imbalance induced mode contaminants. Excitation amplitudes and phase slowness of borehole guided modes are presented for the most common excitation regimes, i.e., monopole, dipole, quadrupole and the hexapole excitation, the latter one showing to have the advantage of a higher formation shear cut-off frequency than the quadrupole mode. Special emphasis will be on the analysis of the dipole excitation and the differences that occur due to variations in tool rigidity, tool diameter and (isotropic) formation properties with the resulting conclusion that a formation flexural mode is not observable as a result of a LWD dipole excitation (this opposed to its dipole wireline counterpart). The guided mode excitation amplitudes are calculated as residues using a Laurent series expansion. This (unconventional) way of calculating the residue has the advantage that it is independent of the pole order, does not require the numerical evaluation of derivatives with respect to the vertical wavenumber and allows for an accurate and efficient FFT implementation.     

Magnetostratigraphy of deep drilling core SG-1 in the western Qaidam Basin (NE Tibetan Plateau) and its tectonic implications

The Qaidam Basin is the largest intermontane basin of the northeastern Tibetan Plateau and contains a continuous Cenozoic sequence of lacustrine sediments. A ~ 1000-m-deep drilling (SG-1) with an average core recovery of ~ 95% was carried out in the depocenter of the Chahansilatu playa (sub-depression) in the western Qaidam Basin, aimed to obtain a high-resolution record of the paleoenvironmental evolution and the erosion history. Stepwise alternating field and thermal demagnetization, together with rock magnetic results, revealed a stable remanent magnetization for most samples, carried by magnetite. The polarity sequence consisted of 16 normal and 15 reverse zones which can be correlated with chrons 1n to 2An of the global geomagnetic polarity time scale. Magnetostratigraphic results date the entire core SG-1 at ~ 2.77 Ma to ~ 0.1 Ma and yielded sediment accumulation rate (SAR) ranging from 26.1 cm/ka to 51.5 cm/ka. Maximum SARs occurred within the intervals of ~ 2.6–2.2 Ma and after ~ 0.8 Ma, indicating two episodes of erosion, which we relate to pulse tectonic uplift of the NE Tibetan Plateau with subsequent global cooling.